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ABSTRACT 

We develop models of void formation starting from a small initial fluctuation at re- 
combination and growing to a realistic present day density profile in agreement with 
observations of voids. The model construction is an extension of previously developed 
algorithms for finding a Lemaitre-Tolman metric that evolves between two profiles of 
either density or velocity specified at two times. Of the 4 profiles of concern — those of 
density and velocity at recombination and at the present day — two can be specified 
and the other two follow from the derived model. 

We find that, in order to reproduce the present-day void density profiles, the initial 
velocity profile is more important than the initial density profile. 

Extrapolation of current CMB observations to the scales relevant to proto-voids 
is very uncertain. Even so, we find that it is very difficult to make both the initial 
density and velocity fluctuation amplitudes small enough, and still obtain a realistic 
void by today. 

Key words: Cosmology: cosmic microwave background; cosmological parameters; 
theory; early Universe; large-scale structure of Universe. 



1 AIM 

Voids are vast regions of the Universe with high negative 
density contrast, which are fundamental parts of the large 
scale structure of the Universe. Although the mean radius 
of voids is 10 Mpc, they contain only few galaxies. Ac- 
cording to the data from the 2 degree Field Galaxy Redshift 
Survey (2dFGRS^) processed by Hoyle and Vogeley (2004), 
about 40% of the volume of the Universe is taken up by 
voids. 

The aim of this paper is to describe the non-linear 
growth of voids out of small initial density and velocity per- 
turbations on a homogeneous background at the moment 
of last scattering. We use the inhomogeneous, spherically 
symmetric dust (Lemaitre-Tolman) model, an exact solu- 
tion of Einstein's equations. This paper also reports on the 
main factors responsible for the formation of voids, and a 
simulation of void evolution is presented. As follows from 
our previous papers (Krasinski and Hellaby 2002, Krasinski 
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and Hellaby 2004), the final state is sensitive not just to 
the amplitude, but also to the exact profile of the initial 
perturbations. So although velocity perturbations of rela- 
tive amplitude (AV/V) around 8 ■ 10~^ were needed in our 
models to reproduce realistic voids, it is still possible that 
other profiles can be found for which a smaller initial veloc- 
ity amplitude will sufiice. 



2 A HISTORICAL OVERVIEW 

The discovery of large-scale cosmic voids became possible 
when astronomers started to measure the distribution of 
galaxies in space. Since William and John Herschel's re- 
searches, i.e. from XIX century, it was known that galaxies 
cluster (for example in the Virgo or Coma cluster). How- 
ever there could be no certainty that these clusters were not 
just caused by the galaxies being projected on the celestial 
sphere. This changed with the publication of the Hubble law 
in 1929. Five years later Tolman and, immediately after. Sen 
studied the stability of the Friedmann models with respect 
to inhomogeneous perturbations, and concluded that they 
are unstable against a rarefaction caused by a negative per- 
turbation of either the initial density or the initial velocity. 
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Consequently there should be condensations as well as un- 
derdense regions in the Universe. 

Nevertheless it was only at the end of 1970s that voids 
were discovered. But the very first sign of their existence had 
been known 20 years earlier. In 1960 Mayall had measured 
the redshift of 50 galaxies in the Coma cluster (Mayall 1960) . 
His survey covered 33 square degrees of the celestial sphere. 
The Coma cluster was also studied by Chincarini and Rood 
in 1975, but it was Gregory and Thompson's (1978) large 
survey that ended with the discovery of a void. Their survey 
covered 260 square degrees of the sky and comprised galaxies 
up to 15th magnitude, which is equivalent to a radial velocity 
of approximately 8000 km/s. 

Also Joeveer and Einasto (1978) observed voids in their 
redshift survey. At the beginning of 1980s the term 'void' 
was first used to call the regions avoided by galaxies (Rood 
1988). 

Next researchers discovered extremely large regions 
avoided by galaxies called supervoids. The classical example 
is the void discovered by Kirshner et al. (1981). From the 
late 1970s, they were observing the galaxies in the Bootes 
and Corona Borealis constellation and discovered an almost 
empty region with the size of 50-100 Mpc, (where h is 
the Hubble parameter in units of 100 km/s/Mpc)). More 
up-to-date measurements determine this size to be approxi- 
mately 60 Mpc and prove that this region is not entirely 
empty. Dey et al. (1990) observed 21 galaxies within this re- 
gion and estimated the mean density contrast to be 

5^ P"--'-P^ e r-0.84. -0.661 (1) 

pmat 

where pmat and pmat are the density and mean density of 
matter in the Universe. 

Other discoveries of voids followed very fast. In 1982, 
after a five year survey, the CfA (Center for Astrophysics) 
galaxy redshift catalogue was finished. Huchra's team mea- 
sured redshifts for over 2400 galaxies with luminosity below 
magnitude 14.5. The amount of data increased in later sur- 
veys. For example, the current Sloan Digital Sky Survey 
(SDSS^) covers one-fourth of the celestial sphere and will 
measure redshifts for over 10** objects. 



3 SIZES AND SHAPES OF THE VOIDS 

The size of a void depends on the luminosity of the galaxies 
that surround it. The research done by Lindner et al. (1995) 
has shown that the mean void size, as estimated by bright 
elliptic galaxies, varies between 13 and 36 Mpc, whereas 
the size evaluated by the measurement of fainter galaxies 
drops to 9 - 25 Mpc. This phenomenon is known as 
the void hierarchy, whose characteristics are listed in Table 
□ (Lindner et al. 1996). 

According to the data from 2dFGRS, the average radii 
of voids in NGP (North Galactic Pole) and in SGP (South 
Galactic Pole) are 14.89 ± 2.67 h'^ Mpc and 15.61 ± 2.48 

Mpc respectively (Hoyle and Vogeley 2004). The max- 
imal radius of a sphere inscribed into a void is 12.09 ± 
1.85 h-^ Mpc in NGP and 12.52 ± 1.99 h'^ Mpc in SGP. 

2 www.sdss.org 



Table 1. The void hierarchy (from Lindner et al.l996). 



Type of object 


Mean 


size 


rich clusters (Abell Catalogue) 


100 h 


-1 Mpc 


poor clusters (Zwicky Catalogue) 


37 h- 


1 Mpc 


bright (M < -20.3) elliptical galaxies 


30 h- 


1 Mpc 


galaxies brighter than M = -20.3 


23 h- 


1 Mpc 


galaxies brighter than M = -19.7 


16 h- 


1 Mpc 


galaxies brighter than M = -18.8 


13 h- 


1 Mpc 



These sizes were estimated using galaxies whose luminosity 
Mum — 5log{h) varied from -18 to -21. 

The computer algorithm employed by Hoyle and Voge- 
ley in their void search used a similar rule to that invented 
by El- Ad and Piran (1997). Firstly, the algorithm decides 
whether a particular galaxy should be assigned to the wall 
(a region of higher density surrounding the void) or to the 
void itself. Then the voids are filled with spheres. 

It is apparent that the shapes of small voids are close to 
spherical, while the largest voids are more irregular. How- 
ever, they can still be divided into smaller spherical voids. 
Sato (1982 and 1984) once hypothesized that spherical voids 
collide with each other as they expand, to produce Zel- 
dovich's "pancakes". The fact that large nonspherical voids 
can be divided into spherical regions may thus indicate that 
they are conglomerates of smaller spherical voids that had 
already collided. Together with another result, that for voids 
the spherical shape is stable (Sato and Maeda 1983), this 
shows that the L-T solution is the right device to model 
voids. 



4 DENSITY CONTRAST 

While considering the above data, the following question 
arises: are these regions really empty? Or are some objects 
just not bright enough to be detected by galaxy surveys? 
This issue was investigated by Thun and co-workers in 1987. 
They examined the distribution of faint, dwarf and irregular 
galaxies from the Nilson UGC catalogue and concluded that 
faint galaxies do not form any separate large-scale structures 
but occur together with bright galaxies. Similar results were 
obtained by Peebles (2001), who analyzed the data from the 
Optical Redshift Survey (ORS). The conclusion is - since 
there are no bright galaxies in voids, there should be no faint 
ones, either. Some observational studies of dwarf galaxies in 
the nearest neighborhood (Pustil'nik et al. 1995; Popescu et 
al. 1997; Kuhn et al. 1997; Gorgin and Geller 1999) found 
galaxies in voids. However, these were only single objects, so 
the deduction of Thun and Peebles was qualitatively right. 
Today voids are defined as regions of low density, or as re- 
gions avoided by certain type of objects, for example bright 
galaxies. According to these definitions, voids do not have 
to be empty. Rojas and co-workers (2003) focused on the 
photometrical properties of galaxies in the voids from the 
SDSS data. They defined the voids Eis regions of density 
contrast S < —0.6. In their study, they measured galaxies 
up to 2 = 0.089. From 155 000 galaxies they extracted a 
sample of 13742, in which 1010 were galaxies in voids. Their 
analysis shows that galaxies in voids are fainter, more blue 
and more compact than galaxies in the walls. 
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Figure 1. The average density contrast for voids in the 2dFGRS, 
in units of void radius (figure taken from Hoyle & Vogeley 2004) . 



Of the galaxies from the 2dFGRS, 5 % are void galax- 
ies. The average density contrast for voids in this survey is 
presented in Fig0 The mean density contrast in the central 
region is 5 = -0.94 ± 0.02 in NGP and 5 = -0.93 ± 0.02 in 
SGP. 

However, data obtained from observations of galaxies 
give information about the luminous matter, while dark 
matter does not have to concentrate in galaxies. Peebles 
(2001), referring to the existing radioastronomical observa- 
tions of 21 cm waves (sensitive to HI clouds) (Weinberg at 
al. 1991; Hoffman, Lu and Salpeter 1992; Szomoru at al. 
1996; Zwaan at al. 1997) and to the results of observations 
of Lyman a absorption line system (Bergeron and Boisse 
1991; Steidel, Dickinson and Persson 1994; Lanzetta et al. 
1995) concluded that not only galaxies, but also gas clouds 
avoid voids. This conclusion can only be justified indirectly 
because there are no direct observations of Lyman a clouds 
in voids. 

The voids that are discovered in surveys are defined 
by galaxies with low redshift. For example, the maximal 
redshift used by Hoyle and Vogeley (2004) to find voids in 
2dFGRS was z — 0.138. For such a small redshift, the Lyman 
a line (A = 1216 A) is in the ultraviolet, which is beyond 
the reach of observations from the ground. The resolution 
of lUE (International Ultraviolet Explorer) satellite was not 
sufficient to detect the Lyman a clouds. 

One must ask whether luminous matter is a good tracer 
of mass, or is there a significant amount of dark matter 
within voids? Is the density contrast for dark matter the 
same as for luminous matter, or do galaxies prefer regions of 
higher density contrast? In the second case, the real density 
contrast would not be as low as the one in Fig. 

Answers to these questions suggested by cold dark mat- 
ter N-body simulations are inconclusive. In these simula- 
tions, the positions of test particles in a chosen part of space 
are traced. Progress in this field depends of the computing 
power. In 2002, Arbabi-Bidgoli and Miiler used 256"^ par- 



ticles of masses IO^^/i'^Mq and 4 • lO^°/i"^M0. They ob- 
tained density contrasts from S ~ —0.6 for the voids with 
the mean diameter 10 - 16 Mpc to 5 ^ —0.8 for the 
voids with mean diameter 36 - 50 Mpc . Benson with 
co-workers (2003) used 512^ particles of masses 10^ h-^MQ. 
Their simulations implied that the density contrast should 
be (5 ~ -0.8, 0.85. The Gottlober team (2003) used 1024^ 
particles of 4- 10^ h~^MQ and obtained the density contrast 
of approximately 5 = —0.9. Some N-body simulations face a 
crucial problem, as they predict that voids should be filled 
with dwarf galaxies. 

In our models, we assume that the real density con- 
trast in voids is the same as the one obtained from galaxy 
observations. 



5 THE LEMAITRE - TOLMAN MODEL 

5.1 Formulae and general properties 

The Lemaitre-Tolman model is a spherically symmetric solu- 
tion of Einstein's equations with a dust source. In comoving 
and synchronous coordinates, the metric is: 



2 J .2 

: c at 



l + 2E(r) 



dr 



i^^(^,r)d!:^^ 



(2) 



where d^^ = de^ + sin 9d(j> , and E{r) is an arbitrary func- 
tion of r. Because of the signature (-I-, — , — , — ), this function 
must obey E{r) > — |. 

The Einstein equations reduce to the following two; 



2M,r 



R 



(3) 



(4) 



where M(r) is another arbitrary function and k — 

When R^r ~ and A/^ 7^ 0, the density becomes infi- 
nite. This happens when shell crossings occur. Equation ||1J 
is similar to its Newtonian counterpart for a spherical dust 
distribution: 

l-r,2 _ f ^ GM 

where TZ, £ and A4 are respectively the radial coordinate, 
the energy of the particles, and the mass within radius TZ (in 
Newtonian mechanics, the cosmological constant is not con- 
sidered). Therefore, M{r)(? /G is the mass inside the shell of 
the radial coordinate r, and E{r)c? is the energy per mass 
unit. 

Equation @ can be solved by simple integration: 



dR 



y 2£ + 2|. + iAi?2 



c[t-tB(r)], 



(5) 



where Ib appears as an integration constant, and is an ar- 
bitrary function of r. This means that the Big Bang is not 
a single event as in the Friedmann models, but occurs at 
different times for different distances from the origin. 

When A = 0, the above equation can be solved in para- 
metric form: 
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For E <0: 



M 



R = - — (l-cosr?), 



77 — sm rj = 



3/2 



M 



i{t-tB{,r)). 



(6) 



Eliminating 77 one can write this as'^ 



cIb 



ct 



M 



(-2^)3/2 



arccos 1 + 2 



, ER 
If 



R 



1 - 

For E = 0: 
lMc\t~tB{r)y 



^ ER 
' M 



(7) 



(8) 



• For £■ > 0: 

R = 

sinh r] — r] = 
or equivalently: 



^(coshry-1), 



{2E) 



3/2 



-c(f-fi3(r)). 



(9) 



Cts 



ct 



M 



(2S)3/2 



1 + 2 



ERV 
M j 



arcosh 1 + 2 



M 



Thus, the evolution of a Lemaitre-Tolman model is deter- 
mined by three arbitrary functions: E{r), M{r) and tB{r). 
The metric and all the formulae are covariant under arbi- 
trary coordinate transformations of the form r = f{r'). Us- 
ing such a transformation, one function can be given a de- 
sired form. Therefore the physical initial data for the evolu- 
tion of the Lemaitre-Tolman model consist of two arbitrary 
functions. 



5.2 The Friedmann limit 

The Lemaitre-Tolman model is a generalization of the Fried- 
mann models and becomes a Friedmann model when the 
following conditions are fulfilled: 



(i) tB = constant. 

This constant is usually set to zero. 

|£|3/2 



(ii) 



M 



constant. 



(11) 



(12) 



In the Friedmann limit, the density distribution is a function 
of the time coordinate only, and is expressed by the formula: 

i^pc = (13) 

The above conditions are invariant under any coordi- 
nate transformation. In the class of coordinates used here, 
one can choose the radial coordinate as: 



^ Equation applies with rj < n, i.e. in the expansion phase. 
We do not consider the recollapse phase in this paper. 



R{r,t)^rS{t), (14) 

where S is the scale factor of the Friedmann models, and 
then 



M(r) = Mor-^ 



E{r) = Ear 



(15) 



The Friedmann limit is an essential element in our ap- 
proach. As mentioned above, our model of void formation 
describes a single void in an expanding Universe. Far away 
from the origin, the density and velocity distributions tend 
to the values which they would have in a Friedmann model. 
The mean sizes of voids presented in the literature are es- 
timated in the Friedmann models (each estimate uses one 
specific model, but for low redshifts the differences between 
different models are negligible). Using 11411 one can identify 
the areal radius of a void R{r^,to) with the mean void radius 
given in the literature. 



6 BACKGROUND MODELS 

The aim of this paper is to describe the formation of voids 
from initial density and velocity perturbations at the time of 
last scattering, and also to check how the evolution of a void 
depends on the background model. This requires knowledge 
of the density and velocity perturbations, and also of the age 
of the Universe at the moment of last scattering. Although 
not all the background models used here are consistent with 
observations, even those excluded by observations will clar- 
ify some mechanisms responsible for void formation. The 
astronomical data put limits on the values of some parame- 
ters. From the observation of the oldest stars the lower limit 
for the age of the Universe is estimated to be approximately 
12 — 14 X 10^ years (Spergel et al. 2003). From the measure- 
ment of the movement of galaxies in clusters and of matter 
in galaxies one can estimate the mean matter density. In 
the critical density units this value is Qmat ~ 0.3. Obser- 
vations of type la supernovae and of the microwave back- 
ground radiation suggest a nonzero cosmological constant, 
of approximate value SIa ~ 0.73 (Bennett et al. 2003). 

In the following, the subscript t indicates a background 
value, indicates a present day value, and u indicates a 
value at last scattering. 



6.1 Background models without cosmological 
constant 

6.1.1 The Einstein-de Sitter Universe 

The Einstein-de Sitter universe is the flat Friedmann model 
filled with matter and without the cosmological constant. 
This model can be obtained from the Lemaitre-Tolman 
model. From I13II and HIH with _E = it follows that the 
density of the homogeneous background is: 

6M 1 

Pb = -ZJPF^ = TTIT^- (16) 



nc'^Rf 



6-KGt^ ■ 



Substituting above the critical density obtained from the 
Hubble constant Ho = 72 km/s/Mpc (Bennett et al. 2003) 
one can get that the present age of the Universe is 9.053 x 10® 
years. The last scattering photons on electrons took place 
when z ~ 1089 (Bennett et al. 2003). It was not a single 
event, but a process extended in time. In what follows, it 
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will be assumed that the last scattering took place when 
z = 1089. At that instant the background density was equal 
to 

{pis)b = {poUl + zf (17) 

so from 1161 we get that the Universe was 252 000 years old. 
The following quantity is a measure of the velocity: 

R,t 



In the flat background this becomes: 



cMi/3 



4 \3 
Set) 



(18) 



(19) 



6.1.2 The hyperbolic background 

The hyperbolic background is the fc < Friedmann model. 
The age of the Universe at the moment of last scattering 
can be calculated using eqs. @. In the Friedmann models 
the factor A-I/E'^^'^ is constant and one can calculate it as- 
suming that the current expansion rate of the homogeneous 
background is given by the Hubble constant: 



Hi 



i?2 



2Ec^ 



2Mc^ 



From this one gets: 
2E 1 2, ^ 



2Ec^ 
Rl 



2/3 



npoC 



(20) 



(21) 



where pcrit ~ SHq / kc* is the density of a fc = Friedmann 
model. Substituting the above result in the first of (jOJ, one 
gets 



Vis 



ar cosh ( 2 



l + n„ 



(22) 



where Qmat = po/ pcrit = KC^po/iH^. From the above for- 
mula, and from 12111 and @ it follows that the age of the 
Universe with fimat = 0.27 when z = 1089 is equal to 477 
000 years. The present age is 11.1 x 10^ years. In the model 
with fimot = 0.391 (this model has a similar present day 
velocity to the model with fimat = 0.27 and SIa = 0.73) one 
can obtain 402 000 years for the instant of last scattering 
and 10.6 x 10^ years for today. The value of the background 
density was calculated from eq. 11711 . 

The background velocity in the hyperbolic model is: 



cAP/3 



2E 2M1/3 



M2/3 



+ 



Rb 



1/3 



(23) 



6.2 The elliptic background 



The elliptic background is the k > Friedmann model. The 
procedure is similar to the one above. To calculate the age 
of the Universe one must use eqs. Q, and the value of rj is 



rjis = arccos 



^niat{l + z) ^ 

The age of the Universe in the model with ^mat = 11 is 76 
000 years at last scattering and 4.6 x 10^ years today. 



6.3 Background models with the cosmological 
constant 

6. 3. 1 The flat background 

When _E = 0, eq. can be solved explicitly: 

R 

dR'- 



V + 3^-^^ 



4 / 3 2M A , , 

3A-~th(^-^— + -I =c(t-t,). 

or equivalently: 



„3 6M . ,2 , /3A 
H = — — smh ct\ — 
A 



It follows that the background density is: 



KC pb = 



A 



sinh 



(24) 



(25) 



(26) 



The age of the Universe can be calculated from the formula 
given above. Substituting the values suggested by astronom- 



ical observations, Qmat = 0.27 and JIa ~ 



= 0.73 we 



obtain that the density is equal to the critical density, the 
Universe is 13.48 x 10^ years old today, while at the moment 
of last scattering the age of the Universe was 484 000 years. 

The background velocity in the flat Universe with the 
cosmological constant is given by the formula: 



(R.i 



cM3 



KC^pb 



— coth'^ 
4 



ct\ 



(27) 



6.3.2 N on- flat background models 

We follow an analogous procedure to the A = case, and 
calculate: 



2E 



-lipcritC (1 — fimat — ^a) 



/ 6 



2/3 



(28) 



where Qa = {A/ kc^) / pcrit = Ac^/3H^. When E ^ 0, eq. 
does not have an analytic solution. Let us denote: 



Rb 



2E 



then eq. can be written: 



a 

J 



dx 



■ at. 



v"+i + i^^' 



(29) 



In the model with flmat = 0.27 and Qa ~ 1.64 one finds 
that the age of the Universe at the moment of decoupling 
was 485 000 years and at present it is 32.46 x 10^ years. The 
background velocity is: 



bb = 



(R.i 



2E 2M1/3 1 Rf 



V A-/2/3 i?6 3 M2/3 



2E „ / I ^V^^ 1 A 6 
^ + 2^-.pbC^j +-A(^— 



2/3 



(30) 



The parameters of the various background models are 
summarized in Table |5| 
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Table 2. The age of the Universe at the present epoch and at 
the moment of last scattering 



angular size f^l 



Model 






Present age of 
the Universe [y] 


Age at last 
scattering [y] 


f^mai — 


1, Ha = 




9.053 X IQS 


252 X 10^ 


f^mat — 


0.27, Ha = 





11.1 X 10'' 


477 X 10^ 


^mat — 


0.391, Oa = 


= 


10.6 X 10^ 


402 X 10^ 


^mat — 


11, Qa = 




4.6 X lO'' 


76 X 10^ 


f^mat — 


0.27, Qa = 


0.73 


13.48 X 10'' 


484 X 10^ 


^mat — 


0.27, Ua = 


1.64 


32.46 X 10^ 


485 X 10^ 



7 PERTURBATIONS AT THE MOMENT OF 
LAST SCATTERING 

The microwave background radiation is a relic from the 
epoch when the Universe was young and hot. When the 
temperature of the Universe dropped below 3000 K, the 
free path of photons became comparable to the Hubble ra- 
dius and radiation stopped interacting with matter. Conse- 
quently, an analysis of this radiation can give us information 
about the state of matter at that moment. The spectrum of 
the cosmic microwave background is the spectrum of the 
black body with the mean temperature of T ~ 2.725 K 
(Mather et al. 1999). More precise observations measure the 
fluctuations with amplitude ST/T ~ 10~^. These are related 
to the density fluctuations at last scattering. Unfortunately, 
this relation is complicated - one must take into account 
several effects. 




Figure 2. 

measured by 



The temperature 
WMAP (based 



fluctuations (AT) 
on data from 



http:/ /lambda. gsfc.nasa.gov/product/map/ I. 

Table 3. Angular diameters of the pre- void region 
Model Ae 

nmat = 1, Ha = 0.246° 
^mat = 0.27, Ha = 0.071° 
n^at = 0.27, f^A = 0.73 0.143° 



7.1 Linear and angular diameters 

The observations of the microwave background radiation 
measure differences in temperature between two points of 
the celestial sphere. The angular resolution of the instru- 
ments is not high. The Wilkinson Microwave Anisotropy 
Probe satellite (WMAP) measures the radiation in five dif- 
ferent ranges and, depending on the range, has the resolution 
ft-om 0.82° to 0.21° (Bennett et al. 2003). The results of the 
temperature measurement are presented in Fig. |21 

The following question arises: what are the linear diam- 
eters of the regions that WMAP can still see as separate? 
And is this resolution sufficient to observe the regions which 
evolved into the currently observed voids? 

The linear diameters can be estimated using the scale 
law of the Friedmann models: 



where 



Ldec 



xir)So 



Xir)Sdec 

where x{^) is the distance in the space t — const: 

r 

X{r) -.^(l/Sit)) J g,.r{t,r')dr'. 



These diameters are related to the angular diameters: 
L = DAe», 



(31) 



(32) 



(33) 



where D is the angular distance; in the Friedmann models 
it is determined by: 



D 



So 
1 + Z 



(34) 




Hd)-- 

and d is given by the formula 



when 
when 
when 



^mat + f^A > 1 
^mat + f^A ~ 1 
^mat + f^A < 1 



d- 



HqSo 



dz' 



^0^(1 + z')'^ + fi™t(i + z')-^ + riA ' 



(35) 



(36) 



where f2c = 1 — f^mat — ^a. Assuming that the moment of 
last scattering was when z = 1089, and that the mean void 
diameter is 25 Mpc, for different models one gets the 
results shown in table |21 

Comparing data from Table |21 with Fig.|21one sees that 
the resolution of the WMAP instruments is not sufficient to 
make direct measurements of the temperature fluctuations 
in the regions that became voids. 



7.2 Initial fluctuations 

As the starting point of further considerations, we need a 
rough estimate of the initial conditions. The procedure for 
estimating the density fluctuations from the observed tem- 
perature fluctuations is complicated. The main contribu- 
tions to the observed temperature perturbation are from 
the intrinsic temperature fluctuations of the emitting fluid, 
the frequency shift as light emerges from a fluctuation in 
the gravitational potential (Sachs- Wolfe effect), and from 
the Doppler effect (motion of the emitter). The moment of 
last scattering also depends on the cosmological background 
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model, but for diflerent models the differences are negligi- 
ble. (For completeness, one should take into account what 
happened to the radiation on the way to the observer: the in- 
tegrated Sachs- Wolfe effect — the effect of gravitational per- 
turbations along the line of sight, the Rees-Sciama effect — 
the influence of changes of the gravitational potential with 
time, the Sunyaev-Zeldovich effect — radiation interacting 
with galaxy clusters, and also the non-linear consequences 
of the photons travelling in an inhomogeneous space.) 

It is found (e.g Padmanabhan 1996) that the intrinsic 
temperature fluctuations obey {AT /T) intrinsic ~ Ap/3p. 
The magnitude of velocity perturbations is usually not 
given*. Fluctuations in the fluid motion of magnitude AV 
(away from uniform expansion) must contribute a Doppler 
term {AT /T)Doppier = AV/c, so observations of AT/T < 
10~^ put an upper limit of 10~^c on AV. On the other 
hand, since the fluctuations are acoustic oscillations, the 
three contributions to {AT/T) observed must be of compa- 
rable magnitude, {AT /T) intrinsic ~ $/3 ~ AV/c. Indeed 
for an oscillating fluid with a relativistic equation of state 
{dp /dp = c?), we must have AV ~ cAp/ip. 

However, the major problem is not with calculations, 
but with the data. The present data are available for scales 
larger than the scale of a single void, and all these calcula- 
tions must rely on extrapolations. 

The WMAP data that was used as the source for Fig.|21 
has such a large scatter for scales around 0.2° that the tem- 
perature fluctuation could be anything from 10"'' to 10~^. 
Using the extrapolation proposed by the WMAP team one 
would get AT/T ^ 2 ■ 10"^ 

To estimate f = AV/Vb — Ab/bb for a present day scale 
of Lo = 12 Mpc, we write L = Lo/{l + z), H = i?o(l + 2)'''", 
Vf, = LH and n = 2/3 for an Einstein-de Sitter model. Thus 
we flnd, for Ha = 72 km/s/Mpc and z = 1089, that 

v^2x 10"''. (37) 



8 FORMATION OF VOIDS 

In this section we will check whether it is possible to evolve 
the voids from a small initial density and velocity pertur- 
bation imposed on a homogeneous background. We will also 
check what is the influence on the structure formation of the 
following factors: 

(i) The shape and the amplitude of the initial perturba- 
tions, 

(ii) The evolution time, 

(iii) The expansion rate, 

(iv) The outflow of mass from central parts of the void. 

The evolution of the void will be calculated with six dif- 
ferent background models. Some of these models, especially 
the elliptic ones with and without cosmological constant, are 
inconsistent with the observations. These are used not in or- 
der to obtain a model of the observed Universe, but to check 
which of the factors listed above are more important in the 

* It is commonly stated by specialists that AV/c ~ 10"^ but we 
have not been guided to a reference that says this or something 
equivalent. 
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process of void formation. For comparing the various mod- 
els, let us assume that the initial conditions are independent 
of the background model. 

To determine the evolution of the Lemaitre-Tolman 
model one needs to know two functions. In this paper, these 
functions will be the initial density and velocity distribu- 
tions. This is not the only method to specify the evolution 
in the Lemaitre-Tolman model. The evolution can be deter- 
mined also by giving the initial and the flnal density pro- 
files (Krasihski and Hellaby 2002) or the initial velocity dis- 
tribution and the final density distribution (Krasihski and 
Hellaby 2004). ^ 



8.1 The algorithm 

The computer algorithm used to calculate the evolution of 
a void was written in Fortran and consisted of the following 
steps. Numerical methods are from Press et al. (1986) and 
Pang (1997). 

(i) The initial time tt was chosen to be the time of last 
scattering. This moment was calculated as described in sec- 
tion 16.3.21 Equation (1291 was integrated using Bode's rule 
with the step of value 10"^ • {Q/ {Kpcr(?i^mat)Y^^ . The ho- 
mogeneous background density and velocity were calculated 
from eqs. Ill7t and 13011 as described in section 16.3.21 

(ii) The initial density and velocity fluctuations, imposed 
on this homogeneous background, were defined by functions 
of radius I, 

5{l) and v{tj, 

as listed in Tables 0] and |S] and the actual density and ve- 
locity followed from 

p^{e■) ^ {pb)^{l + m) and b{l)^{bb)^{l + v{l)), 

with p measured in units of lO'^MQ/kpc'' and b in kpc"''^^. 
The parameter I is defined as the areal radius at the moment 
of last scattering, measured in kiloparsecs, and is also used 
for the radial coordinate, i.e. 

r = l = Ri/d = R{r,ti)/d 

where d — 1 kpc. 

(iii) Then the mass inside the shell of radius Ri , measured 
in kiloparsecs, was calculated by integrating eq. Q: 



M{£) - M(0) = 




Since the density distribution has no singularities or zeros 
over extended regions, it was assumed that Imin ~ and 
i\f = at ^ = 0. 

The integration was done using Bode's rule, with step size 
2.5 pc. 

(iv) E was calculated from Ri, V = {R^t)i, M and a cho- 
sen A value, using eq. ((4j. 

(v) Then ts was calculated from eq. Q using Simpson 
integration, with step size 10"^^. 

^ The numerical examples in those papers used the present day 
for the final time. 
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(vi) Once M, E and Ib are known, the state of the 
Lemaitre-Tolman model can be calculated for any instant. 
Solving eq. Q with the second-order Runge-Kutta method 
for R{t,£) along each constant £ worldline, we calculated the 
value of R{t,£) and R^t{t,£) up to the present epoch. The 
time step was 5 x 10^ years. 

(vii) The density p{t, £) was then found from ^ using the 
five-points differentiating formula. The adjusted differences 
between adjacent worldlines, used in estimating derivatives, 
was 10 pc. 

(viii) The density contrasts presented in section IF!^ were 
estimated from eq. 11391 . 



8.2 The Void Models 

8.2.1 Initial perturbations of homogeneity 

Because of lack of precise observational data, it is not pos- 
sible to calculate the exact profile of the initial density and 
velocity perturbations. From the measurements of the mi- 
crowave background radiation one can estimate only the am- 
plitudes of these profiles. It can be intuitively expected that 
the region which in the future would become a void should 
have, at the initial instant, a minimum of density and a 
maximum of velocity at the center.^ 

The chosen initial density and velocity distributions ful- 
filling the above conditions are presented in Fig. |3 These 
profiles conform to the amplitudes estimated in sec. 17.21 
They were defined by the functions presented in Tables 2] 
and El 

From the current observations, only the average density 
contrast is known. For the purpose of comparing our results 
with the observational data (Fig. , the results shown in 
Fig.2](and also the figures showing the final density contrast 
in what follows) do not present the real density contrast, but 



the average one, i.e.: 



where p is the present background density: 



<p>-- 



(39) 



(40) 



Unfortunately there are no astronomical data for the 
current velocity distribution in the void. It is more practi- 
cal to measure the expansion rate by the equivalent of the 
Hubble parameter (Rt/R). The results are presented in Fig. 

m 

Fig- m shows the re-distribution of mass resulting from 
void formation. Curve O is the function M{R) at the initial 
time t = tis (i.e. M(r) versus R{r,tis)), and the other curves 
show the calculated M{R) at the present time, t = to, for 
different background models. 

Curve 6 in Fig.|l]is truncated. At the cutoff point a shell 



^ But see the papers by Mustapha and Hellaby (2001) and by 
Krasinski and Hellaby (2002) — maxima and minima can be re- 
versed during evolution, and it is not at all necessary that a void 
begins with a minimum of density at the center. 

The form of these initial density and velocity distributions is 
a result of consecutive adjustments made in order to test the 
influence on void formation of the various factors mentioned at 
the beginning of this section. 
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Figure 3. The initial density (upper) and velocity (lower) per- 
turbations for data discussed in section [8.2.1l 



crossing occurs - the inner shells catch up with the outer 
shells. This results in a singularity that probably does not 
occur in the real Universe. Before it happens, the gradient of 
pressure would become significant and the Lemaitre-Tolman 
model would become invalid. 

From the figures presented above some initial conclu- 
sions can be made. One can say that the current depth of 
the void depends of the amount of mass re-distribution, i.e. 
the mass outflow. This outflow depends on the expansion 
rate and on the age of the void. The major influence on the 
structure formation is from the shell expansion. 

In models 4 and 5 the expansion rate is bigger than in 
models 1 — 3 and the final density contrast in models 4 and 
5 is much more negative than in models 1 — 3, even though 
the age of the Universe in these models is much lower. 

Unfortunately, models 1 — 5 cannot recover the ob- 
served density contrast of today's voids. The proper depth 
can be obtained only in model 6, where the age of the Uni- 
verse is significantly larger. 

The age of the void (in limits estimated by the various 
astronomical observations) is of lesser importance, compared 
to the expansion rate. In models 2 and 3 the expansion of 
the shells is similar and the final density contrast in model 
3 is lower due to the time available for evolution being by 
2 X 10^ years longer. 

We now check how big the influence of the initial shape 
of the density and velocity perturbations is. 
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Table 4. The initial Density Perturbations used in the runs. All the values in the table are dimensionless, and the distance parameter is 
the areal radius in kiloparsecs I = i?,;/lkpc. Note that the output figures depend on the initial perturbation in both density and velocity. 

Section Model Density perturbation Parameters Graph Output 
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Fig.iroi 



The explicit formulae for the profiles are presented in Tables 
Hand El 

The final results are not very different from the one 
shown in Fig. 2] This proves that the velocity distribution 
in the formation of voids is very significant, while the density 
distribution is of lesser importance. 



8.2.4 Amplitude 

In this subsection, the amplitude of the initial fiuctuations 
is increased compared to the one used in subsection 18.2.11 
and is 4 ■ IQ-'^. The profiles of the initial perturbations are 
presented in Tables 3 and 4, and the final results in Fig. [71 
The increased amplitude of the initial perturbations re- 
sults in a void with a higher negative density contrast. To 
obtain a density contrast near S ~ —0.9 we needed to in- 
crease the amplitude of the initial density profile more than 
70 times, and the amplitude of the velocity profile 20 times, 
compared to the values estimated from CMB fiuctuations. 
Even so, the value S = —0.94 of the density contrast in the 
void was not reached, except in the two non-realistic back- 
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8.2.2 Homogeneous velocity profile 

For this case, the initial density profile 5 is as in Fig.|31 while 
the initial velocity profile is f = 0. The explicit formulae for 
these profiles are presented in Tables |11 and |3 

The results shown in Fig. |S1 seem to be surprising. The 
mass re-distribution is almost the same in all the models 
(except model 6), the diameters are similar, but the current 
density profiles are different. In this case, a second factor re- 
sponsible for the void formation is seen, which is the faster 
expansion rate. The faster expansion of the void compared 
to the homogeneous background causes that the difference 
between the density in the central regions of the void and the 
density in the background increases with time, even when 
the mass of the shell inside the region of R{t, r) is not chang- 
ing very much. Consequently, in the models with higher ex- 
pansion rate the density contrast is most negative. 

8.2.3 Homogeneous density profile 

In contrast to the above, we now set the initial density profile 
to (5 = 0, while the initial velocity profile, v, is as in Fig. Rl 
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Table 5. The Initial Velocity Perturbations used in the runs. All the values in the table are dimensionless, and the distance parameter 
is the areal radius in kiloparsecs t = Ri/l kpc. Note that the output figures depend on the initial perturbation in both density and 
velocity. 
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Figure 4. The current density contrast and the Hubble parame- 
ter for data discussed in section r8.2.1l in six different background 
models: 1 - n^at = 0.27, Oa = 0; 2 - n,nat = 0.39, = 0; 
3 - Omat = 0.27, Ha = 0.73; 4 - Q^at = 1, = 0; 5 - 
Umat = 11, ^A = 0; 6 - Umat = 0.27, f^A = 1-64. * - shell 
crossing. 




Figure 5. The mass re-distribution for data discussed in section 
18.2.11 Ro is the smallest R values at which the density takes the 
background value. Mq is the mass inside the shell of areal radius 
Ro- — the initial condition; other labels as in Fig. HI 
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Figure 6. The current density contrast and the Hubble param- 
eter for the flat initial velocity profile in different background 
models. Labels as in Fig. HI 



ground models. In model 4 the minimum value is —0.908 
and in model 3 it is —0.873. Unfortunately, increasing the 
amplitude leads to a shell crossing singularity in some mod- 
els. 



8.2.5 Observation and a model - a cross check 

In the previous sections we had problems generating voids 
from small initial density and velocity fluctuations. The only 
alternative was to use a background model with an ex- 
tremely large age of the Universe (inconsistent with limits 
estimated by the various astronomical observations). 

In this subsection, we will try to choose the initial pro- 
files that lead to the best fit with observational data. We 
focus only on one background model, preferred by the astro- 
nomical observations, which is Q.mat = 0.27 and 0,^ = 0.73. 
It should be noted that the difference between void evolution 
in this model and in the Einstein-de Sitter model is not big, 
so the problem of void formation is in the initial conditions 
rather than in the background model. 

So far none of the results obtained have recovered the 
observational data, presented in Fig. Q This figure presents 
the density profile of the void with very smooth edges. Un- 
fortunately this profile does not reach to regions where the 
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Figure 7. The current density contrast and the Hubble param- 
eter for the initial data with higher amplitude of density and ve- 
locity perturbation, in six different backgrounds models. Labels 
as in Fig. |1] 



density is higher than the background density, and where the 
superclusters of galaxies would be found. The mean value of 
the density contrast inside the void is 5 ~ —0.94. The den- 
sity contrasts of the voids obtained so far have been too 
shallow or had steep edges. 

The initial fluctuations are presented in Fig.|H| The pro- 
files are presented in Tables 01 and |5] The results are shown 
in Fig. l^and llUI Model 2 has both proper density contrast 
and smooth edges. The conclusion from numerical experi- 
ments with different shapes of the initial profiles is that a 
model of a void consistent with observational data (with the 
value of density contrast less then 5 — —0.94, smooth edges 
and high density in the surrounding regions) is very hard 
to obtain within the Lemaitre-Tolman model, without the 
occurrence of a shell crossing singularity. The final state of 
model 2 was very close to this singularity and in model 4 a 
shell crossing occurred. 

The main factor responsible for void formation is the 
velocity perturbation, with an amplitude of ~ 8 ■ 10~^ , near 
the centre, and dropping below zero in the outer regions 
(model 3 did not fulfill this condition). The density fluctua- 
tion is of lesser importance. The models 1 — 3 had the same 
initial density fiuctuations and in model 4 the amplitude 
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Figure 8. The initial density (upper) and velocity (lower) per- 
turbations for the results presented in Figs. 151 and IIUI 



was two times greater. In spite of these differences, the final 
results differ only in shape, and not in the depth of the final 
density contrast. 

In simulations with higher values of background den- 
sity the initial velocity and density perturbation needed to 
obtain similar final results are smaller (see subsection 18 . 2 .41 . 



9 EVOLUTION 

Let us take a closer look at the evolution of model 2. 

Fig. ll2l shows the density distribution in eight different 
moments of time. In this section the figures do not present 
the density contrast, but the real density distribution calcu- 
lated from eq. Q. 

Figs. ITHl - 1151 show respectively the functions M{R), 
tsiR) and E{R), where R is the areal radius at the initial 
instant (curves O) and at the final instant (curves 1). The 
pictures demonstrate the evolution of the structure: in the 
expanding void mass moves outwards. 

As one can see from Fig. 1141 the difference in the value of 
the Bang time function between outer and central regions of 
the void is of the order of hundreds of years. This is negligible 
compared to the age of the Universe which at the moment 
of last scattering was of the order of 10^ years. 
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Figure 9. The current density contrast for the models discussed 
in section [8.2.5l 
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Figure 10. The Hubble parameter for the models discussed in 
section 
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Figure 11. Comparison of the curves from Fig.|5]{2,4) with the 
observed density contrast from Hoyle and Vogeley (2004) (SGP 
and NGP). 




Figure 12. The evolution of the density distribution in a void 
for the age of the Universe: 10 (No. 1), 100 (No. 2) million years 
and 1 (No. 3), 5 (No. 4), 10 (No. 5) billion years after the Big 
Bang, and the current (No. 6) profile of density distribution for 
the model discussed in section |5] -Ro is the point at which the 
density takes the background value for the first time. 




Figure 13. The mass re-distribution for data discussed in section 
1^ Ro is the point at which the density takes the background value 
for the first time. Mq is the mass inside the shell of areal radius 
-Ro- Curve O is the function M(R) with R and Rq taken at ti, 
curve 1 is the function M{R) with R and Rq taken at time t2. 



10 LIMITATIONS ON THE INITIAL 
CONDITIONS 

An exact reconstruction of the initial conditions that held 
at the time of decoupling in the region that would become 
a void is not possible due to the lack of precise data. Nowa- 
days the only way to reconstruct the initial conditions is the 
theoretical approach. The linear theory is often used for this 
purpose because of its simplicity. 



10.1 Initial conditions in the linear theory 

In linear theory the shape of the density perturbation is 
constant in time. Only the amplitude is changing, according 
to the formula: 



5{t) = DSo, 



(41) 
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Figure 14. The Bang time tg as a function of R{ti , r) (curve 
O) and of R{t2, r) (curve 1) for the model discussed in section 1^ 




Figure 15. The energy function E as a, function of R{ti,r) 
(curve O) and of R{t2,r) (curve 1) for the model discussed in 
section 1^ 



where So is the current value and D is the linear growth 
factor. In the Einstein-de Sitter Universe it is equal to; 

1 



D = 



1 + z' 

and in the general case: 

a 

D{a) = [ f{a')da', 



2a/(a) 



where 
/(«) = 



^l + nM(i-l)+nA(a2-l)' 



while a = 1/(1 + z). For large redshifts in any model: 
Dia) ^ 



1 + z 



(42) 



(43) 



(44) 



(45) 



In the linear regime in the spherically symmetric case, 
the peculiar velocity of the density perturbation is given by 
the formula (Peebles 1993, pages: 115-116): 



(46) 



where H is the Hubble constant, S is the scale factor and / 
is the velocity factor. For large redshifts in any cosmological 
model this asymptotically becomes unity. Since: 



rr _ S,t 



R.t = rS.i 



(47) 



the velocity fluctuations are determined by the formula: 
^ = -^5 (48) 

That implies that to calculate the evolution of the void one 
needs the initial density and velocity perturbations with the 
amplitudes: 



: 9 ■ 10" 



u^3 10~ 



(49) 



These results are more than 10 times smaller than in the 
Lemaitre-Tolman model, where the amplitude needed is 8 • 
10"^ 

Since the current evolution is not linear, one cannot 
use the linear theory to predict the initial conditions in the 
young Universe. However the initial conditions could be ob- 
tained with the help of the Lemaitre-Tolman model, since it 
is an exact solution of the Einstein filed equations. There still 
remains the question about its accuracy, since the Lemaitre- 
Tolman Universe is a very simple model of a void. One can 
try to answer this question by a careful look at the restric- 
tions of the Lemaitre-Tolman model. 



10.2 Main limitations of the Lemaitre-Tolman 
model 

The Lemaitre-Tolman model is spherically symmetric, so it 
cannot take rotation into account. The density and velocity 
distributions can only depend on the radial coordinate. How- 
ever, astronomical observations show that for voids these 
conditions are fulfilled with a satisfactory accuracy. 

The next limitation is the dust energy-momentum ten- 
sor, which means that the pressure of matter and of radi- 
ation are neglected. At the present epoch this is correct, 
but one can ask about the error at the moment of the last 
scattering implied by this assumption. For the purpose of 
estimating that error, a general Friedmann - Lemaitre - 
Robertson - Walker (FLRW) model will be used. 

(i) The pressure: 
The equation of state for the perfect gas is: 



p= , (50) 

where fcs is the Boltzmann constant and mn is the proton 
mass. Because of lack of data about the nature of dark mat- 
ter, let us assume that the mean molecule weight /x = 1. In 
FLRW models, T = ro(l -I- 2), where To is the present value 
of the background temperature. Substituting the numerical 
values we obtain that the ratio of the pressure to the density 
of matter is: 



^ 2.75 • 10-^0. 



(51) 



(ii) The radiation energy density: 



The expression for the density of the radiation energy is: 



€rad = aT =aTo{l + z) , 



(52) 
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where a = — , and cr is Stefan-Boltzmann constant. The 
density of matter in the FLRW models is; 

p^Po{l + zf. (53) 

If we assume that that the present value of density is equal 
to the critical density in the flat FLRW model, then the 
ratio of the radiation energy density to the matter energy 
density will be: 



^rad 



aTl 



+ ^0.054 



(54) 



From the above one can see that the gas pressure is 
negligible. Ignoring the radiation one makes a small error 
which at first seems to be of little importance. However, from 
the above simulations with different background models, it 
can be seen that if the initial energy density is higher then 
the final contrast of matter density is deeper. Unfortunately, 
this leads to a shorter age of the Universe. To obtain the 
proper age of the Universe, a higher value of the cosmological 
constant is needed. 

Radiation avoids this problem because during the evo- 
lution it becomes less and less significant and the age of the 
Universe does not change so much. Is it possible then that a 
model of void formation that includes radiation could pre- 
dict the voids observed today starting from smaller initial 
perturbations? 



10.2.1 Radiation 

The Einstein field equations for the spherically symmetric 
perfect fluid distribution can be reduced to the two following 
equations (Lemaitre 1933): 

K,R^R,rPC = 2Af ,., 



KR^R,tp = -2M,t, 



(55) 
(56) 



where p is the energy density, while p is the pressure. M is 
defined by the formula: 

2M(r,t) = R{r,t)R/{r,t) - 2E(r)R(r,t) - iAii^(r, t).(57) 

As in the case without radiation, M(? /G is equal to the 
mass inside the shell of radial coordinate r. In this case, the 
mass is not constant in time and in the expanding universe 
it decreases. 

From the equations of motion T°''^-^ = we obtain 

P,r = 0, 
P,e = 0, 
P,^ = 0, 



P.t + 



P 



R Rn 



0. 



(58) 
(59) 



Equations 1)58^ require that pressure can only be a function 
of time. If inhomogeneous radiation should be included in 
the model, the metric form should be changed — the goo 
component should be a function of the radial coordinate. 
This would require finding a new exact solution of the Ein- 
stein equations. 

However, it is instructive to know what changes are 
caused by homogeneous radiation. The time dependence is 
the same as in the Friedmann models: 



10 12 14 16 18 20 22 24 

R [Mpc] 




10 12 14 16 18 20 22 24 

R [Mpc] 



Figure 16. The final density contrast and the Hubble parameter 
in the model with (+) and without (x) the radiation. 



^rad 



tde 



i/3 



Prad — g^rad- 

For the purpose of comparing with the previous results, the 
background model is the flat Friedmann model with A = 0. 
The initial density and velocity distributions are presented 
in Table H and m 

The computer algorithm implemented to do the calcula- 
tion was similar to the one used in models without radiation. 
The only difference is that instead of solving one equation 
(eq. ||1J), the second-order Runge-Kutta method was used 
to solve simultaneously eqs. (1561 and 15711 . 

Fig. ll6l shows the comparison between the models with 
and without radiation. The final density contrast in the 
model with radiation is a little lower inside the void, but 
higher at the edge compared to the model without radia- 
tion. 

As one can see, homogeneous radiation is not of great 
importance in the formation of voids. The difference from 
the model without radiation is not big because, in spite of 
higher energy density, as one can see from eq. (I5(ill . the mass 
of the shell is decreasing. The decrease of the mass implies, 
via eq. 15711 with A = 0, that the velocity of the shell is 
also decreasing with time. Consequently the evolution of 
the structure slows down and the final results do not show 
significant differences compared to the model without radi- 
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ation. These results suggest that the error estimated with 
eq. 1541 is indeed of little importance and even the model 
with inhomogeneous radiation should not lead to significant 
differences in the final results. 



11 CONCLUSION 

The main aim of this paper was to produce a non-linear 
model of void formation, starting from small initial density 
and velocity fluctuations that existed at last scattering, and 
to investigate what factors are necessary to reproduce cur- 
rent observations. 

• In the numerical experiments that we carried out, the 
perturbations of that were needed to form a realistic present 
day void had to have a density amplitude of 5p/ p ~ 5 • 10"'^ 
and a velocity amplitude 5V/V ~ 8 ■ 10"'^ (in the model with 
Qmat = 0.27 and f^A = 0.73). 

• It was found that density perturbations are of lesser sig- 
nificance than velocity perturbations in the process of void 
formation. 

• There was no significant difference between the evolu- 
tion of the void in the model with and without the cosmo- 
logical constant. In this latter case the amplitude of the ini- 
tial velocity fluctuation needed for this purpose was a little 
smaller. 

• The existence of voids is closely related to the existence 
of regions of higher density surrounding the voids. In our 
simulations, there were problems obtaining reasonable pro- 
files for these high density regions, since shell crossing singu- 
larities tended to occur. The Lemaitre-Tolman model breaks 
down at shell crossings, so to trace the further evolution of 
the void we had to focus on the central regions of the void. 
(In reality, as density increases, a gradient of pressure would 
appear, which cannot be described in the Lemaitre-Tolman 
model.) Since superclusters are observed on the edges of the 
voids, we interpret this singularity as an indicator of the 
presence of a supercluster. 

Thus our numerical experiments were not entirely suc- 
cessful at generating a void consistent with observational 
data, out of perturbations at last scattering that would be 
as small as indicated by current structure formation theo- 
ries. Either the initial perturbations had to be larger than 
the conventional values, 5p/ p ~ 10~^ ~ SV/c, or else the 
present-day void was too shallow. This discrepancy calls for 
an explanation. In attempts to explain it, several hypothe- 
ses might be considered. We list a set of such hypotheses in 
the order of decreasing probability of being correct; this is 
of course our subjective evaluation. 

(i) The experiments, here and in Paper II (Krasinski and 
Hellaby 2004), showed that the final state is sensitive not 
just to the amplitude of the initial velocity perturbation, 
but also to its profile (the density perturbation is less sig- 
nificant). We may not have identified the profile that gives 
the best consistency with observations. Work on this will be 
continued. We tried one approach already, but the results 
were not encouraging. We took the density distribution in 
the void at present as given, and the density distribution 
at recombination homogeneous. We calculated the implied 
velocity distribution at recombination, which was 40 times 



too large (see Sec. 18. 2. St . Then we numerically decreased 
the values of that velocity by the factor of 40, took it as 
part of the input data, and calculated the implied density 
distribution at recombination. It was again too large, so we 
decreased it numerically by the appropriate factor, took it 
as input, and calculated the implied velocity field at recom- 
bination. This iteration quickly converged to stable values of 
the amplitudes: 6 • 10~^ for density and 7 • 10"'^ for velocity. 
Clearly, this does not solve our main problem. 

(ii) Voids may not be as empty as they appear. It is possi- 
ble they contain a significant amount of unobserved matter, 
such as gas or other baryonic "dark matter" . A present day 
density contrast of smaller absolute value is easier to pro- 
duce with small initial fiuctuations. 

(iii) Matter may have more components than comoving 
dust, and these other components may be dynamically sig- 
nificant. Including other components into a fully non-linear 
description, and allowing for the possibility that the various 
kinds of matter do not comove, requires finding solutions of 
Einstein's equations with two or more independent streams 
of matter as a source, which is an extremely challenging 
problem. Some possible matter components are: 

(a) Cold dark matter, i.e. some form of non-baryonic 
matter that decouples from normal matter very early on 
and starts forming structures before last scattering. The 
amplitude of such fluctuations could be larger than what 
CMB measurements allow for baryonic matter. Whilst 
there is evidence for some non-luminous matter from 
galaxy rotation curves, galactic interactions and gravita- 
tional lensing, it is less than the required "concordance" 
value, and of unknown composition. The difficulty with 
this hypothesis is that cold dark matter is not based on 
any confirmed physical theory and has yet to be detected. 

(b) A radiation component that is still significant just 
after last scattering. 

(iv) Observations give the Galaxy distribution in redshift 
space, and void sizes are deduced by using a Hubble law 
based on a homogeneous Friedmann model. Since the den- 
sity distribution thus obtained is clearly not homogeneous, 
the use of a Friedmann model is really not correct. This may 
introduce quite significant errors in the density and velocity 
profiles of voids. The well known "finger of God" effect is an 
example of a large discrepancy between true positions, and 
the Friedmann-based mapping from redshift space to physi- 
cal distances. Consistency with Einstein's equations requires 
use of an inhomogeneous cosmological model to correctly 
map redshift space into physical distances and thus deter- 
mine how velocities vary with distance. At present, neither 
the density distribution nor the velocity field of galaxies is 
known with confidence.** Had we known both these fields 
with a reasonable precision, we might use the L-T model to 
calculate the fluctuations in density and velocity at recom- 
bination. 

(v) Reliable measurements of temperature fluctuations 
of the cosmic microwave background (CMB) radiation are 
available only for angular scales larger than 0.5° (see Fig,|2J, 



° The velocities are measured, but then, without the inconsistent 
assumption of homogeneity, we do not know to which points in 
the curved manifold of the Universe they should be attached. 
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with the scatter becoming large at 0.2°. On the other hand, 
as shown in Paper I (Krasihski and Hellaby 2004) and in 
sec. 17.11 the angular diameters, in the CMB sky, of struc- 
tures the right size to evolve into voids, are less than 0.25°. 
Though current theory suggests the amplitude decreases at 
smaller scales, this has yet to be confirmed by improved 
CMB temperature measurements. 

(vi) Solid observational data give the temperature fluc- 
tuation of the cosmic microwave background radiation at 
present, AT/T, separately for different modes of pertur- 
bation. However, there are several factors contributing to 
AT/T (see section [7.21 . which may partially cancel each 
other. We were interested in the magnitude of the fluid ve- 
locity at last scattering, as an initial condition for our mod- 
els, but we were not able to locate a formula relating the 
fluid velocity to the observed AT/T. Thus our estimate of 
AV may be inaccurate. We ourselves do not plan to enter 
this field, however we would welcome an explicit calculation 
with a reliable result. 

(vii) The Universe may be much older than currently be- 
lieved, making the time available for void formation much 
longer. This would require adjustment of the matter content 
of the universe, of the value of the cosmological constant or 
of the value of the Hubble constant, each of which may affect 
the structure-formation timescale. 

(viii) General Relativity (GR) may not be the right the- 
ory of the evolution of the Universe. There is certainly no 
lack of proposed modifications or alternatives. While some 
would gladly accept this conclusion, and no-one can claim 
with certainty that GR will survive all future scrutiny, we 
wish to stress that GR has a much stronger experimental 
basis than any cosmic structure formation theory. 
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